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ABSTRACT 

A simulation of the thermonuclear explosion of a Chandrasekhar-mass C+O white dwarf, the most popular 
scenario of a type la supernova (SN la), is presented. The underlying modeling is pursued in a self-consistent 
way, treating the combustion wave as a turbulent deflagration using well tested methods developed for labo- 
ratory combustion and based on the concept of 'large eddy simulations' (LES). Such consistency requires to 
capture the onset of the turbulent cascade on resolved scales. This is achieved by computing the dynamical 
evolution on a 1024 3 moving grid, which resulted in the best-resolved three-dimensional SN la simulation car- 
ried out thus far, reaching the limits of what can be done on present supercomputers. Consequently, the model 
has no free parameters other than the initial conditions at the onset of the explosion, and therefore it has con- 
siderable predictive power. Our main objective is to determine to which extent such a simulation can account 
for the observations of normal SNe la. Guided by previous simulations with less resolution and a less sophis- 
ticated flame model, initial conditions were chosen that yield a reasonably strong explosion and a sufficient 
amount of radioactive nickel for a bright display. We show that observables are indeed matched to a reasonable 
degree. In particular, good agreement is found with the light curves of normal SNe la. Moreover, the model 
reproduces the general features of the abundance stratification as inferred from the analysis of spectra. This 
indicates that it captures the main features of the explosion mechanism of SNe la. However, we also show that 
even a seemingly best-choice pure deflagration model has shortcomings that indicate the need for a different 
mode of nuclear burning at late times, perhaps the transition to a detonation at low density. 

Subject headings: Stars: supernovae: general - Hydrodynamics - Instabilities - Turbulence - Methods: nu- 
merical 



> 

(N 
O 

o 
o 



1. INTRODUCTION 

Understanding the physics of Type la supernova (SN la) 
explosions is a necessary premise for using them as cosmic 
distance indicators with confidence. Consequently, the de- 
velopment of three-dimensional explosion models has been a 
major challenge in astrophysics over the past decade. Only re- 
cently, however, a self-consistent three-dimensional treatment 
of the turbulent thermonuclear flame in a white dwarf (WD) 
star avoiding free parameters became available. The ultimate 
goal of these efforts is to describe the full observational set of 
SNe la and even to predict their observed properties, such as 
their luminosities, light curves, and spectra, from first princi- 
ples, with sufficient accuracy to guide the empirical calibra- 
tion of their peak luminosity. 

While this stage of sophistication is clearly not yet reached 
with current simulations, at least one (perhaps in some 
cases the dominant) part of the mechanism can already 
be treated in an non-parameterized approach. In the first 
stage of the supernova explosion the flame propagates 
subsonically and is accelerated by turbulence creating a 
problem of turbulent combustion physics. The knowl- 
edge on this subject acquired for laboratory combustion 
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in engineering sciences can be transferred to astrophysics, 
facilitating self-consistent simulations with particular em- 
phasis on th e correct representation of the fl ame/turbulence 
interacti on (iNiemever & Hillebrandtl 119951: IReinecke et al.1 
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Rtipke & Hillebrandt 
I2006bt ISchmidt & Niemeyerl 12006: 
ISchmidt et al.l l2006bl) . In fact, most of the properties of 
turbulent chemical flames are recovered in SNe la, and 
the largely different relevant length scales, the simpler 
reaction kinetics, and the absence of complicated boundaries, 
make computations of thermonuclear supernovae perhaps 
easier than designing a new car engine, provided the initial 
conditions of the explosion are known. 

Recent results from s everal three-dimensi o nal simulations 
(IReinecke et al.l l2002bl: iGamezo et al.l 120031: iTravaglio et al 
| 2004t IRopke & Hillebrandtl l2005at iBravo & Garcfa-Sen: 
boOfjt IRopke et al l l2006bt ISchmidt & Niemeyerl l2006l) ~iri- 
dicate that a turbulent thermonuclear deflagration — a sub- 
sonic flame propagation — alone can lead to powerful ex- 
plosions of the WD star. In combination with the 
elimination of free parameters in the model the ques- 
tion arises to what extent the results of such simula- 
tions match observational constraints. This question mo- 
tivated several previous studies which addressed initial pa- 
rameters of the exploding WD (IRopke & Hillebrandtl l2004t 
Trava glio et al.1 120051: IRopke et al.l l2006al), the flame igni- 
tion configuration (IReinecke et al.l2002bURopke et all2 006b ; 
Schmi dt & Niemever 2006), s ymmetry con straints in the sim- 
ulation setup (IRopke & Hi llebrandt 2005a), a nd the p articular 
choice of the turbulence model (Schmidt et al. 2006b), as well 
as the impact of some initi al condiditions on observational re- 
sults dBlinnikov et al.ll2006l) . 
One question left out until now, however, was whether the 
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results would change significantly with higher resolution. The 
flame/turbulence interaction determines the propagation of 
the burning front and thus the results of the model with respect 
to nucleosynthesis and energetics. Given the expected strong 
turbulence effects and the associated vast range in relevant 
length scales this endeavor may seem hopeless. However, as 
discussed below, with the numerical modeling strategy pur- 
sued here, it turns out to be sufficient to resolve the onset 
of turbulence on the computational grid in order to recover 
a definitive result. The necessary resolution can be reached 
with current computational resources and marks the aim of 
the present study. A direct comparison with observational 
data is possible and is indeed pursued in this publication for 
one particular simulation in the framework of the pure de- 
flagration model. With what is so far the best resolved three- 
dimensional simulation of the deflagration scenario for SNe la 
we assess its consistency and its compatibility with observa- 
tions. As we show in this paper, such a state-of-the-art defla- 
gration model matches observational data surprisingly well, 
although not all details are reproduced perfectly. As is ar- 
gued below, the simulation reaches a state where the physical 
effects are resolved so that its predictions are expected to be 
robust and to reflect generic properties of the deflagration sce- 
nario. 

Of course, a single simulation cannot provide the basis for 
general conclusions on the validity of the deflagration model. 
Here, we present a simulation performed within the range of 
realistic initial conditions (as far as known) but explicitly set 
up the initial flame configuration in a favorable way with re- 
gard to energy and 56 Ni production. Combining this single 
realization of the deflagration model with the exploration of 
the parameter space carried out previously (and partially with 
simpler modeling approaches and less resolved simulations), 
successes and limitations of the deflagration model become 
visible and will be discussed accordingly (see Sect.[7j. In par- 
ticular, some persisting shortcoming of the deflagration model 
can be taken as an indication that the deflagration stage alone 
cannot explain all aspects of SNe la. In some cases models be- 
yond the de flagration regime are supported by obse rvational 
results (e.g. iMazzali et all 120071; IWang et al.l 120071) . but un- 
certain or even hypothetical physical mechanisms generally 
introduce free parameters into the models and make a face-to- 
face comparison with observations more difficult than in the 
case of the pure deflagration model. 

2. ASTROPHYSICAL MODEL 

Recent observations (e.g. Mazza lTet alj|2007l) support the 
scenario that the majority of SNe la originate from a 
carbo n-oxygen WD d isrupted by a thermonuclear explo- 
sion dHovle & Fowled Tl 9601) when it approaches the Chan- 
drasekhar mass limit by mass accretion from a binary 
companion which may pe rhaps be a non-degenerate star 
(Ruiz-Lap uente et all 12004); for alternative progenitor sce- 
narios see e.g. iHillebrandt & Niemeverl |2000b . Caused by 
temperature fluctuations inside the convective core of the 
electron-degenerate WD, small regions near its center are 
thought to commence a thermonuclear runaway with rapidly 
increasing nuclear energy release from burning the carbon- 
oxygen fuel to heavier elements. This process marks the birth 
of the thermonuclear flame. The exact number and distribu- 
tion of ignition kernels is not well determined yet, implying 
that the largest uncertainty lies in the initial parameters of the 
explosion. The virtue of three-dimensional explosion mod- 
els is that they directly relate the ignition conditions to the 



outcome of the explosion without introducing additional pa- 
rameters. Thus a comparison with observations sheds light on 
two aspects of the supernova problem — the progenitor evolu- 
tion giving rise to these initial conditions and the explosion 
mechanism itself. 

Here, we are concerned with the latter, for which several 
possibilities have been suggested. Hydrodynamically, two 
fundamentally distinct modes of flame propagation are admis- 
sible: a subsonic deflagration and a supersonic detonation. A 
prompt detonation would incinerate the entire star at the high 
initial densities, leading to conditions where nucle ar reactions 
produ ce almost exclusively iron group elements (lArnett et alJ 
11971 . These are found to be abundant in the ejecta, predom- 
inantly produced as radioactive 56 Ni. In its decay to 56 Co 
(and later 56 Fe) gamma-rays and positrons are emitted which 
down-scatter to optical wavelengths in the SN ejecta and give 
rise to the observed optical display. At later times, in the so- 
called "nebular" spectra, this iron can be directly observed. 
However, failure to produce the intermediate mass elements 
(e.g. silicon, calcium, and sulfur) observed with apprecia- 
ble abjmd^ncesjnflie early spectra of SNe la (for a review 
see Filippenko 1997) rules out the prompt detonation mecha- 
nism. In contrast, a subsonic deflagration leaves enough time 
for the stellar materia l to expand before i t is reached by the 
thermonuclear flame (Nomot o et al.l 1 19761) . Therefore, some 
burning takes place at densities low enough such that mostly 
intermediate mass elements are synthesized, while most 56 Ni 
is produced during the early deflagration phase. 

A laminar deflagration wave, mediated by microphysical 
heat transport, would propagate far too slowly to explode the 
WD. However, as burning ignites near the WD's center and 
propagates outward, a stratification of low-density ashes be- 
low denser fuel emerges which is buoyancy-unstable in the 
gravitational field of the star. This configuration eventually 
leads to the formation of ascending burning bubbles. Their in- 
terfaces are subject to shear. Since typical Reynolds numbers 
of the flow are of the order of 10 14 , strong turbulent eddies 
must form which decay to smaller scales in a turbulent cas- 
cade. This has tremendous impact on the flame propagation. 
The flame surface is corrugated by turbulent velocity fluc- 
tuations and its area increases dramatically, which enhances 
the burning rate. This acceleration of the flame speed and 
the fuel consumption rate enable explosions of the WD (e.g. 
iReinecke et alll2002bl) . 

The energy input into turbulent eddies is due to the 
Rayleigh-Taylor instability on the largest scales (^lOOkrn). 
However, the flame interacts with turbulent velocity fluctu- 
ations down to much smaller scales. As demonstrated in 
Sect. below, these velocity fluctuations originate from an 
inertia-driven turbulent cascade. They do not obey a sim- 
ple buoyancy-dominated scaling and therefore require special 
modeling effort. 

3. NUMERICAL IMPLEMENTATION 

Numerical schemes to model the propagation of a deflagra- 
tion front in a SN la have been developed over the past years 
and have reached a high level of sophistication. The major 
challenge in the numerical implementation of the concept out- 
lined in Sect. |2]is posed by the vast range of relevant scales. 
It spans about 1 1 orders of magnitude, from the radius of the 
WD (~ 2000km and expanding in the explosion) down to the 
Kolmogorov scale (less than a millimeter) where turbulent en- 
ergy is dissipated. The flame interacts with turbulence down 
to the Gibson scale at which the laminar flame propagation 
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speed becomes comparable to the turbulent velocity fluctua- 
tions. For most parts of the supernova explosion this scale is 
considerably larger than the flame width (< 1 mm). There- 
fore, the interaction of the flame with turbulence is purely 
kinematic and burning takes place i n the so-call ed flamelet 
regime of turbulent combustion (e.g. Peters 2000). Here, the 
flame propagation speed completely decouples from the mi- 
crophysics of the burning and is entirely determined by tur- 
bulence effects. This also implies that the flame propagation 
does not scale with the laminar burning speed. 

The challenge of scales can be tackled in an elegant way by 
employing a large eddy simulation (LES) strategy. This ap- 
proach resolves only the largest scales of the turbulent cascade 
and models the unresolved scales, so that it becomes possible 
to accommodate the entire WD on the computational domain. 
When seen from the scale of the WD, the flame appears as 
a sharp discontinuity separating the fuel from the ashes. A 
suitable numerical scheme to track t he propagation of such 
an interface is the level s et method jO sher & Sethianl [l988t 
ISmiljanovski et all Il997t iReinecke et alJ Il999bl) . where the 
flame front is associated with the zero level set of a signed 
distance function G set up to be negative in the fuel regions 
and positive in the ashes. The propagation velocity of this in- 
terface is given by the effective turbulent burning speed at the 
grid scale which in the flamelet regime is completely deter- 
mined by turbulent velocity fluctuations on that scale. These, 
in turn, are derived from a subgri d-scale turbulence model (for 
details see Sch midt et alJl2006alfbT) which follows the proper- 
ties of turbulence on unresolved scales based on filtering the 
velocity fluctuations on the resolved scales. This way, a local- 
ized treatment of the turbulent properties of the flow becomes 
feasible and the flame/turbulence interaction is modeled self- 
consistently. 

The hydrodynamics is described by the reactive Euler 
equations (with gravity as external force) and numerically — 
discretized on a Cartesi an computational g rid — treated in a 
finite volume approach dFryxell et al.ll 1989b . Of course, the 
tremendous expansion of the WD during the explosion has to 
be taken into account. This was achieved by employing two 
nested moving grids — a fine uniform inner gri d tracking the 
flame and a coarse o uter grid capturing the WD ('Ropke l2005t 
iRopke et alf2006bl) . 

With this implementation we performed what is thus far 
the largest simulation of a thermonuclear supernova explo- 
sion. It was set up on more than a billion computational 
grid cells (a 1024 3 cells grid). This allowed for a very fine 
resolution (down to ~ 870 m) of the central part of the WD 
at the beginning of the explosion facilitating an exception- 
ally fine-structured initial flame configuration that was com- 
posed of 1600 spherical kernels of radius 2.6km distributed 
inside a sphere extending 180km around the center of the 
WD (see Fig. [T] upper left panel). Although the details of 
the ignition process are yet unknown, analytical and numer- 
ical st udies support such a multi-spot ignition configu ration 
(e.g. lGarcia-Senz & Woosleylll995HWoosley et alJl2004 and 
the flame ignition kernel s are expected to be o f the order of 
a kilometer in diameter (Iapi chino et al.l 120061) . The config- 
uration chosen here is certainly one that favors burning near 
the center and (as m ulti-spot ignition scenarios do in general, 
Ropke et al. 2006b) enhances the burning and energy release 
over that found in centrally or sparsely ignited scenarios. 

The nuclear reac tions were treated in th e approximative ap- 
proach outlined by IReinecke e t al. (2002a). Only five species 
were taken into account: the fuel consisted of a mixture of car- 



bon and oxygen in equal amounts by mass, intermediate mass 
elements were represented by magnesium, and a mixture of 
56 Ni and alpha-particles represented nuclear statistical equi- 
librium (NSE) material. Depending on the fuel density ahead 
of the flame, p u , the material was converted to NSE (if p u > 
5.25 x 10 7 gcm~ 3 ), in which case the nickel/alpha ratio was 
adjusted according to temperature and density, to intermediate 
mass elements (if 5.25 x 10 7 gcm~ 3 > p u > 5.0 x 10 5 gcm" 3 ), 
or left unburned (for lower fuel densities). The threshold 
for burning to intermediate mass elements was chosen lowe r 
than in previous simulations. As argued by Schmidt (2007), 
the level-set description of the flame propagation may be ex- 
tended beyond the onset of distributed burning at which tur- 
bulent eddies penetrate the internal structure (the "thin reac- 
tion zones regime"). This commences as soon as the flame 
width becomes larger than the Gibson scale — an effect that 
is expected to occur at fuel densities ^10 7 gcm" 3 . However, 
at densities below ~40 6 gcm" 3 nuclear reactions will be slow 
compared to the expansion time scales and, moreover, tur- 
bulence will interact with the burning zone inside the flame 
structure giving rise to the completely distributed ("broken 
reaction zones") regime. Therefore, the burning description 
adopted here may be too optimistic at low densities and over- 
predict the total amount of intermediate mass elements syn- 
thesized in the explosion. On the other hand, at the late stages 
where low fuel densities are reached, the amplitude of tur- 
bulent velocity fluctuations decreases rapidly. Therefore not 
much material is burned here, such that the fuel consumption 
at very low densities becomes insignificant and the potential 
error introduced by uncertainties is expected to be small. 

4. EVOLUTION OF THE EXPLOSION PHASE 

The evolution of the explosion stage proceeds in agreement 
with the expectations outlined in Sect. Q] After ignition in a 
multitude of only partially overlapping flame kernels (see up- 
per left panel of Fig. [TJ, the ash region grows due to burning 
and starts floating towards the surface. The individual kernels 
deform from their initially spherical shape into "mushroom 
cap"-like structures — typical for the non-linear stage of the 
Rayleigh-Taylor instability — and grow. This causes them to 
merge, so that a connected flame structure forms. This struc- 
ture evolves dominated by buoyancy instability on large scales 
and driven by the turbulent cascade on smaller lengths. The 
continued development of substructure and the merger of fea- 
tures create a deflagration structure with a complex pattern 
(see upper right panel of Fig. |T}. Burning and flotation drive 
the flame towards the surface of the WD. The fuel density 
drops as the flame moves outwards and due to radial strati- 
fication and the overall expansion of the WD caused by the 
energy deposit from nuclear burning. Once the fuel density 
falls below the threshold for the production of intermediate 
mass elements, nuclear burning ceases. This occurs first at the 
leading features of the flame and subsequently in more cen- 
tral flame regions. Finally, no burning takes place anymore. 
At this point the outer ash features have reached the surface 
layers of the ejecta (lower left panel of Fig. [T]). The follow- 
ing seconds in the evolution are characterized by the hydro- 
dynamical relaxation towards homologous expansion. This 
stage i s reached to g ood approximation at about 10 s after ig- 
nition (Ropke 2005). 

This three-dimensional evolution leads to a remnant of the 
explosion with characteristic properties. The density structure 
has patterns from unstable and turbulent flame propagation 
imprinted on it (cf. lower right panel of Fig.[T|) and ash regions 
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FIG. 1. — Evolution of the thermonuclear supernova explosion simulation. The zero level set associated with the thermonuclear flame is shown as a blue 
isosurface and the extent of the WD is indicated by the volume rendering of the density. The upper left panel shows the initial set up and the close-up illustrates 
the chosen flame ignition configuration. The subsequent two panels illustrate the propagation of the turbulent flame through the WD and the density structure of 
the remnant is shown in the lower right panel. 



extend to the outermost layers of the expanding cloud of gas. 

The chemical composition of the ejecta of the explosion as 
followed in the hydrodynamical simulation is shown in Fig. [2] 
The outer parts consist mainly of unburned carbon and oxy- 
gen, but some "bubbles" of intermediate mass elements are 
also present near the surface (Fig. [2] left). The central part of 
the ejecta is dominated by iron group elements (Fig.[2]right). 
Inside the ash bubbles some mixing of different species is 
found, but still regions dominated by either iron group ele- 
ments or intermediate mass elements can be identified. The 
large-scale bubbles of ash are separated by a few thin but 
dense fingers of unburned material so that the spherically av- 
eraged composition of the ejecta appears strongly mixed. 

5. CONSISTENCY OF THE SUBGRID-SCALE 
MODELING 

The subgrid-scale model employed in the simulation cou- 
ples the dynamics at unresolved scales to the numerically 
computed flow. The model is based on the balance equa- 
tion for the subgrid-scale turbulent energy, where the en- 
ergy production rate follows from a similarity assumption 
linking the smallest r esolved to the largest unresolved scales 
dSchmidt et al.ll2006ah . Although the calculation of the rate of 



turbulence energy production is localized and, consequently, 
does not rely on stationarity or large-scale homogeneity, 
asymptotic isotropy towards smaller length scales is required 
for consistency. In the case of incompressible turbulence, 
asymptotic isotropy becomes manifest in the Kolmogorov 
A;" 5 / 3 scaling for sufficiently large wavenumbers k. This scal- 
ing has been foun d in simulations of Ra yleigh-Taylor unsta- 
ble stratifications (Cabot & Cook 2006) as we ll as in buoy- 
ancy u nstable flames in degenerate WD matter IZingale et all 

dlop). 

Fig. |3] shows compensated energy spectrum functions, 
k^'^Eik), derived from our simulation for several instants of 
time ranging from 0.5 to 1.0s after ignition. For turbulent 
velocity fields obeying Kolmogorov scaling, the compensated 
spectrum functions are expected to be constant. Wavenum- 
bers are normalized by 7r/(512Ao), where Ao is the cell size 
in the uniform part of the grid. Since the grid is inhomo- 
geneous and does not satisfy periodic boundary conditions, 
the velocity fields obtained from the simulation were zero- 
padded in the outer regions with Gaussian windowing func- 
tions in order to compute Fourier transforms. This is why 
the compensated spectrum functions are plotted for normal- 
ized wavenumbers greater than 16 only. At all times consid- 
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FIG. 2. — Composition of the remnant after the explosion (f = 10s). The products of mass fractions of three of the species followed in the hydrodynamical 
explosion simulation with density are volume rendered in different colors. The right panel shows a cut-away to illustrate the composition of the central parts. 
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FIG. 3. — Compensated energy spectrum functions derived from our simu- 
lation at different instants in time. 

ered here, it can be seen that the spectrum functions are well 
approximated by the Kolmogorov power law for intermedi- 
ate wave numbers, thus indicating the presence of an iner- 
tial subrange. Therefore, the simulation is as well resolved 
as the modeling approach requires. Higher resolution is not 
expected to change the characteristics of the turbulent flame 
propagation. 

6. RESULTS AND COMPARISON WITH 
OBSERVATIONS 

The present deflagration simulation of a thermonuclear su- 
pernova makes distinct predictions for observables that can 
be tested against real SNe la. With certain observables like 
polarimetry spectra and spectra of the nebular phase it may 
be possible to check details of the results of pure deflagration 
models against observations and potentially assess its validity. 
This certainly warrants further attention and will be addressed 
in future work. Here, we are concerned with some more basic 
global properties which test whether our particular simulation 
falls into the range of observations. 

The asymptotic kinetic energy amounted to 8.1 x 10 5() erg. 
In the hydrodynamical simulation, 0.606M Q of iron group 
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FIG. 4. — Spherically averaged composition resulting from the hydrody- 
namical explosion simulation (solid lines) compared to the findings of the 
abundance tomography of SN 2002bo (dotted lines). Iron group element 
abundances are shown in red, intermediate mass elements in green and un- 
burned material in blue. 

elements and 0.43 IMq of intermediate mass elements were 
produced, leaving behind 0.368 Mq of unburned carbon and 
oxygen. 

The composition of the ejecta can be directly compared 
with that of observed supernov ae obtained from abundance to- 
mography (Steh le et alj |2005) using a spherical average over 
the entire ejecta of the model in homologous expansion. Such 
a comparison with the abundance tomography of SN 2002bo 
is shown in Fig. |U where the velocity is used as radial co- 
ordinate. Two cautionary remarks are necessary here. First, 
SN 2002bo was a much more energetic and luminous event 
than that expected to result from the simulation. Second, as 
pointed out in Sect. @l there exist features with high concen- 
trations of particular species (see Fig.|2| large enough to dom- 
inate spectral observations. The procedure of spherically av- 
eraging the mass fractions washes out these structures and ar- 
tificially increases the mixing of species. Nonetheless, the 
chemical composition of the ejecta of the simulated explo- 
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sion and the abundance tomography of SN 2002bo compare 
surprisingly well. In both cases, the central parts are clearly 
dominated by iron group elements and the overall distribution 
of the main species in our simulation generally agrees with 
that obtained from the abundance tomography. There is good 
agreement up to velocities of ~10,000kms -1 , if we keep in 
mind that SN 2002bo had a 56 Ni content of ~ O.5M , which 
is significantly larger than in the simulation. The slope of 
the decline of the iron group abundance towards larger radii 
is also well-reproduced. These elements are mixed out to the 
same radii in the simulation and in the observation. Moreover, 
both agree in intermediate mass elements being more abun- 
dant than unburned material below ~10,000kms _1 and mixed 
compositions over a wide range of velocities. The abundance 
of unburned material in the simulation is low in this velocity 
ra nge, in particu l ar nea r the center. As previously pointed out 
bv lKozma et alj (120051) . the presence of large amounts of un- 
burned fuel at low velocities gives rise to a pronounced oxy- 
gen line in the late-time ("nebular") spectra which is inconsis- 
tent with observations. Owing to the more complete burning 
of the center with the chosen flame ignition condition, this 
problem is alleviated in the current simulation. The abun- 
dance of fuel material near the center is lower by a factor of 
10 to 50 than the iron group mass fraction here. 

Still, the abundance of carbon and oxygen seems rather high 
in comparison with the observation. This may in part be an 
artifact of the spherical averaging procedure for determining 
the composition of the ejecta in the simulation, overestimat- 
ing the contribution of the dense but narrow fuel regions (see 
Fig. 0. As the volume of ejecta that influences the spectra 
taken form a certain line of sight varies with time, a careful 
analysis of this effect would require the derivation of spectra 
from our simulation in a three-dimensional approach. 

At velocities above ~40,000kms _I , however, there is a sig- 
nificant discrepancy between the abundance tomography data 
of SN 2002bo and the simulated composition. While in the 
former the intermediate mass elements dominate out to high 
velocities, the ejecta in the simulation quickly become dom- 
inated by unburned material. This discrepancy, however, is 
only found in the outer 0.25 Mq. 

In order to derive light curves from the simulation, more 
detailed information about the chemical composition of the 
ejecta is necessary than the data obtained directly from 
the hydrodynamical explosion simulation. This informa- 
tion can be derived in a nucleosynthetic postprocessing step 
(Travaglio et al.l 120041) on the basis of 53 3 tracer particles 
that have been advected in the flow. The mass of the star 
was equally distributed between those tracers and they add 
a "Lagrangian component" to the otherwise Eulerian ap- 
proach. Recording representative temperature and density 
profiles they facilitate the a posteriori reconstruction of the 
nuclear reactions with a large reaction network (which would 
be impractical to employ concurrently with the explosion sim- 
ulation). This way, the chemical composition of the ejecta as 
needed for the radiative transfer is determined. In particular, 
it is necessary to separate the 56 Ni, which by radioactive de- 
cay provides the energy for the optical event, from the stable 
iron group isotopes that do not contribute to the luminosity 
but are important for absorption. The postprocessing yields 
0.56M Q of iron group elements (0.05M Q less than the es- 
timate from the hydrodynamical simulation) and 0.33 M of 
56 Ni. This indicates that the simulated supernova would range 
among the lower end of observed normal SNe la. For events 
like SN 1991M, SN 1993L, SN 2004eo, and SN 200 lei simi- 
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FIG. 5. — Bolometric light curve derived for our model (black curves; solid 
is the "UVOIR-bolometric" light curve and the complete bolometric light 
curve is dott ed.) The blue dotted curves correspond to observed bolometric 
light curves (Stritzinger et al. 2006). 

lar masses of S6 N i and total iron group elements were inferred 
bv lMazzali elail (12007b . 

On the basis of the ejecta composition as derived from 
the postprocessing step and spherically averaged over the 
entire star, bolometric light curves (true bolometric as well 
as "UVOIR-bolometric") have been derived with the ra- 
diative transpor t code Stella, as described in detail by 
iBlinnikove t al. (2006). The result is shown in Fig. [5] and 
compared t o the bolometric light c urves of standard SNe la 
as given by Stritzin ger et al.l (|2006). Evidently, the synthetic 
light curve derived from the simulation is a normal SN la light 
curve both in terms of the peak luminosity and the shape. It is 
slightly di mmer than SN 1992A, a s expected from the lower 
56 Ni mass (IStritzinger et all 120061 report O.4M of 56 Ni for 
SN 1992 A), but has a very similar rise and decline rate. 

7. LIMITS OF THE PURE DEFLAGRATION MODEL 

The limited success of the presented simulation in re- 
producing the observables of SNe la raises the question of 
whether it is possible to assess generic limits of the pure de- 
flagration model on this basis. 

As a realization in a single simulation does not 
cover the full parameter space, this has to be dis- 
cussed in combination w it h simu lat ions performed pre- 
viousl y. iReinecke et ail d2002bl): IRopke & H illebrandt 
J2004I); iTravaglio et al l "(120051) : IRopke et al.l (l2006allbh . and 
ISchmidt &^ Niemeverl (120061) explored the impact of various 
initial conditions on the explosion process. These simulations 
were not performed with the same resolution as the one pre- 
sented here and are partially based on the simpler subgrid- 
scale t urbulence model proposed by Niemev er & Hilleb randt 
d 19951) . Still these results can guide an assessment of the de- 
flagration model since the changes between the two subgrid- 
scale models are moderate (at most 30% in the iron group 
elemen t production and the energy release; Schm idt et al.1 
2006b]). 

The two main shortcomings identified in Sect. [6] were (i) the 
overall low explosion strength resulting in a low 56 Ni produc- 
tion and a low explosion energy and (ii) the underproduction 
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of intermediate mass elements in the outer layers of the ejecta. 

Problem (i) has been noted in all deflagration simulations 
based on a consistent approach of modeling the flame prop- 
agation. There are, however, several initial parameters that 
may remedy the problem. These have been tested in previous 
simulations. 

Clearly, a major parameter is the geometrical configura- 
tion of the fla me ignition. With simple central ignitions 
("C3"- setups of iReinecke et al.ll2002bt IRopke & Hillebrandtl 
l2005al) the asymptotic kinetic energy of the ejecta falls around 
4.5 x 10 erg and the masses of produced iron group elements 
hardly exceed 0.5M© with ~ O.3M of it being 56 Ni. This 
would account for the dimmest examples among the "nor- 
mal" SNe la. Ignition in multiple spots around the center 
evidently increases the explosi on strength. Fo r a mo del ig- 
nited in 30 bubbles per octant, Travagli o et al.l d2004) report 
an asymptotic kinetic energy of about 6.5 x 10 5 Q erg and a 
56 Ni production of O.418M . IRopke et al.l d2006bl) explored 
the impact of multi-spot ignitions in a systematic approach 
and found an optimum of such configurations resulting in an 
asymptotic kinetic energy of 6.8 x 10 50 erg and 0.67 M© of 
iron group elements. The latter value may increase in model- 
ing the ignition as a st ochastic process in time to ~ 0.75 M© 
dSchmidt & Niem ever 2006). Given these values, the choice 
of the ignition configuration in the present model was moder- 
ately optimistic producing 0.606M© of iron group elements, 
about 0.33 M© of which is 56 Ni. A large ratio of produced 
stable to radioactive iron group elements (here it amounts to 
0.59) seems to be a general feature of turbulent deflagration 
models ignited in a multitude of ignition kernels since most 
of the iron group synthesis proceeds here at the early stages 
where high densities prevail such that electron capture reac- 
tions are significant. 

Other parameters that have an impact on the 56 Ni produc- 
tion and the energy release are the central density of the WD 
at ignition and its metallicity. A higher central density natu- 
rally increases the amount of material available that poten- 
tially can burn to iron group elements. Moreover, the in- 
creased gravitational acceleration will enhance the buoyancy- 
induced turbulence and accelerate the flame. On the other 
hand, with increased densities electron captures become more 
important and will favor the production of stable iron group 
isotopes over 56 Ni. The effect of varying the central den- 
sity has been addressed by Rtipke et al. (20063). In a density 
range where electron captures are still moderate, increasing 
the central density by a factor of 2.6 changes the 56 Ni yield 
only by less than 10% although the total iron group element 
synthesis was enhanced by ^30%. Therefore, although en- 
ergetic results and enhanced production of the iron group are 
expected with higher central densities of the exploding WD 
star, as electron capture reactions become important the frac- 
tion of 56 Ni in the NSE material and thus the brightness of 
the corresponding events will decrease. An overproduction of 
stable iron group nuclei is anticipated here. Turning this ar- 
gument around, a way to obtain higher 56 Ni masses (although 
the total mass of iron group and intermediate mass elements 
would not increase) would be to consider a WD with lower 
central density ignited in a multi-spot scenario. This issue 
will be addressed in a forthcoming study, but as the results of 
Ropke ^tail (|2006a) indicate the effect is not expected to be 
sufficient to account for the brighter observed SN la events. 

The metallicity of the WD material will have little effect 
on the iron group production and the energy release, but evi- 
dently affects the ratio of stable to radioactive iron group ele- 



of neutron-rich stable nuclei over 56 Ni ( 


Timmes et all 12003: 


Travaglio et al. 2005; Ropke et al. 2006b 


) and makes the cor- 



responding event dimmer. The 56 Ni masses decrease linearly 
with metallicity and the values given above were derived un- 
der the assumption of solar metallicity. Thus, metallicity ef- 
fects have no potential of increasing the 56 Ni yields of the 
deflagration model significantly. 

Issue (ii) is equivalent with the incomplete burning of the 
outer layers of the WD star. Here, the problem of modeling 
the flame propagation at conditions where the flamelet pic- 
ture breaks down becomes important. It has beens shown that 
a flame acceleration in this phase would significantly enhance 
the pro duction of intermediate mass elem ents and t he ene rgy 
release IRopke & Hillebrandl (l2005bl) . but lSchmidtl(l2007l) ar- 
gues that such an acceleration does not occur and that a con- 
tinuation of modeling the flame propagation as in the flamelet 
regime is possible. This approach was followed here, too. 
Therefore, within the pure turbulent deflagration regime it 
seems unrealistic to enhance burning at low densities beyond 
what was achieved in the current simulation. 

We thus conclude that the initial parameters of the model 
are unlikely to leave enough room to adjust the simulations 
in order to account for the more powerful SN la events and 
to reproduce the chemical composition in the outer layers of 
the ejecta. As the model of flame propagation itself leaves 
no free parameters, only such initial parameters of the WD 
and the flame ignition come into consideration here. These 
have been tested and adjusting them in a favorable way allows 
to enter the range of the weaker "normal" SNe la in terms 
of 56 Ni production, explosion energy, and composition of the 
inner part of the ejecta. The more p owerful observed eve nts, 
producing close to 1M© of 56 Ni (e.g. lMazzali et al.ll2007l) are 
clearly out of reach for the pure deflagration scenario of SNe 
la. 

8. CONCLUSIONS 

We presented a three-dimensional deflagration simulation 
of a thermonuclear supernova explosion with the astrophysi- 
cal modeling and the numerical implementation designed in a 
way to allow for self-consistent simulation. The model con- 
tains no tunable parameters except for the flame ignition con- 
dition, which is physically not well determined yet, and the 
initial conditions of the WD which may vary from event to 
event in nature. An additional uncertainty of the model is due 
to burning at the lowest fuel densities which needs further at- 
tention, but is expected to lead to only minor changes in the 
conclusions about the gross properties of the simulation which 
we present here. Clearly, the most significant parameter is the 
flame ignition configuration deciding even about the success 
of t he deflagration burning to explode the white dwarf sta r 
(e.g. lCalder et al.l2004tlPlewa et al.ll2004tlRopke et~aT1l2007l) . 
Here, it was set up in a favorable way in order to burn the fuel 
material at the center of the star as completely as possible and 
to maximize the energy and 56 Ni production. 

The goal of the simulation was to reach a regime in 
which the flame/turbulence interaction modeling is fully self- 
consistent which requires to resolve the onset of the turbulent 
cascade on the computational grid. This was achieved by per- 
forming the simulation with unprecedented resolution. Such 
a consistent implementation of the burning model allows the 
direct confrontation with observables. The comparison with 
bolometric light curves and the chemical composition of the 
ejecta derived from observed spectra showed that the simu- 
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lated explosion indeed reproduces the characteristics of the 
observed fainter examples of normal SNe la. Gross prop- 
erties of the simulations such as energy the energy release 
and the 56 Ni production were also consistent with expecta- 
tions from SN la observations. One has to note, however, that 
even choosing favorable ignition conditions (a slight improve- 
ment is possible with a time-dependent stocha stic distribution 
of ignition kernels, Schmidt & Niemeyer 2006), the explosion 
energy and the brightness of the model event fall towards the 
lowest values observed for normal SNe la. 

Thus, upon comparison with previous results from lower 
resolved simulations exploring the range of the initial param- 
eters, we conclude that within the astrophysical model em- 
ployed here, it seems impossible to reproduce the more ener- 
getic events of the observed SN la sample. Also, a smooth 
shell of intermediate mass elements surrounding the iron-rich 
core of the ejecta cloud and the anti-correlation between the 
iron-group and the i ntermediate mass elements in the ejecta 
(Mazz ali et alj|20 07) are difficult to explain within this frame- 
work. Therefore an extension of the scenario seems necessary. 
The shortcomings of the model in reproducing the chemical 
composition of the ejecta in the outer parts while providing 
good agreement with the observations in central regions may 
taken as an indication that the deflagration phase is only the 
first phase of the explosion mechanism. While the reason- 
able agreement with observables of weaker SN la events indi- 
cates that it may be the dominant part of the explosion mech- 
anism here, the discrepancy is obvious for the brighter parts 
of the SN la sample. Possible extensions of the model include 
an improved description of burning in the distributed burn- 
ing regime at low fuel densities as well as a transition to a 
supersonic detonation mod e of flame propagation at t he on- 
set of this regime. Recently, Ropke & Niemeyer (2007) mod- 



eled delayed detonation scenarios in three-dimensions under 
the assumption that the deflagration-to-detonation transition 
occurs as the flame enters the distributed burning regime. De- 
pending on the configuration of the deflagration flame igni- 
tion, a wide range of 56 Ni masses was recovered which would 
reproduce observations of normal to bright SNe la. In this 
context, the fainter examples correspond to events with a pro- 
nounced deflagration phase. For these, most of the energy and 
56 Ni production originates from that phase and the character 
of the events is similar to pure deflagration models. Thus, 
our results may be interpreted as an indication that the fainter 
sub-sample of normal SNe la expose the hallmark of the de- 
flagration phase. Since the strength of deflagration phase de- 
termined the overall characteristics of the delayed detonation 
model implemented this way, a careful and self-consistent 
treatment of the burning physics in the deflagration phases is 
a premise for credible modeling here as well. 
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